Zeros of the Potts Model Partition Function on Sierpinski Graphs 
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We calculate zeros of the q-state Potts model partition function on m'th-iterate Sierpinski graphs, 
S m , in the variable q and in a temperature-like variable, y. We infer some asymptotic properties of 
the loci of zeros in the limit m — > oo and relate these to thermodynamic properties of the q-state 
Potts ferromagnet and antiferromagnet on the Sierpinski gasket fractal, Soo- 
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I. INTRODUCTION 

Studies of iterated transformations and fractals have 
produced many interesting results in physics and mathe- 
matics (l| . There have been several analyses of spin mod- 
els and critical phenomena on fractals; some early works 
are Q-0- Analyses of zeros of partition functions of spin 
models in temperature-dependent Boltzmann variables 
(Fisher zeros) and field-dependent variables (Yang-Lee 
zeros) have been performed, often for the Ising model, 
and, in some cases, also the Potts model 

Q-df. A dif- 
ferent but related area of interest is vertex colorings of 
one-parameter families of graphs G m that lead to fractals 
in the limit m — > oo, denoted Goo = linim-^oo G m . For an 
arbitrary graph G, the number of ways of assigning q col- 
ors to the vertices of G subject to the condition that no 
two adjacent vertices have the same color is enumerated 
by the chromatic polynomial, P{G, q). 

Here we report calculations of chromatic polynomials 
and their zeros for m'th iterates of Sierpinski graphs, 
denoted S rn , and discuss their asymptotic behavior in the 
limit m — ¥ oo, where 5*00 is the Sierpinski gasket fractal. 
More generally, we also analyze the zeros of the q-state 
Potts model partition function in the complex q plane at 
finite temperature for both the antiferromagnetic (AFM) 
and ferromagnetic (FM) cases. We also give some results 
for zeros in the plane of a temperature-like Boltzmann 
variable (i.e., Fisher zeros). 

We first recall some relevant background. A graph 
G = (V, E) is denned by its set V of vertices (sites) and 
its set E of edges (bonds). We denote n = n(G) = \V\ 
and e(G) = \E\ as the number of vertices and edges of 
G. At temperature T, the (zero-field) Potts model par- 



tition function is given by Z = ^2 
Hamiltonian uM 
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with the 



(1.1) 



where i and j label vertices of G, o~i are classical spin 
variables on these vertices, taking values in the set I q = 
{1, ...,(?}, /3 = (fc^T) -1 , and (ij) label pairs of adjacent 



vertices, or equivalently, edges in E. It will be convenient 
to introduce the notation 

K = /3J , h = /3H , y = e K , v = y-l. (1.2) 

Thus, the physical ranges of v are v > for the Potts fer- 
romagnet and — 1 < v < for the Potts antiferromagnet. 

A spanning subgraph of G is G' = (V, E 1 ) with E' C E. 
The number of connected components of G' is denoted 
k(G'). The Potts partition function can be expressed as 
the sum over spanning subgraphs [16| 



Z(G,q,v) 
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(1.3) 



As is evident from this expression, the partition function 
of the q-state Potts model is a polynomial in q and v 
with positive integer coefficients. For the Potts FM, this 
expression allows one to generalize q from Z + to K+ while 
maintaining a positive Gibbs measure. Eq. (jl.3l) also 
allows one to generalize q and v to complex numbers, 
as is necessary to analyze zeros of Z(G,q,v). For real 
v, the zeros of Z{G, q, v) in the q plane are invariant 
under complex conjugation, and for real q, the zeros of 
Z(G,q,v) in the v or y plane are invariant under complex 
conjugation. 

For the Potts AFM, T -> means K -> -oo and 
thus v — > — 1. As is clear from Eq. (|l.ip . the only spin 
configurations that contribute to Z(G,q,v) in the limit 
v — > —1 are those for which the spins on adjacent vertices 
are different. This proves that Z(G, q, — 1) = P(G,q), 
where P(G, q) is the chromatic polynomial, which, for 
q £ Z + , counts the number of ways of assigning q colors 
to the vertices of G subject to the condition that no two 
adjacent vertices have the same color (called proper q- 
colorings of G). The minimum integer q that allows a 
proper q-coloring of G is the chromatic number, x(G). 

Besides its intrinsic interest in mathematical graph 
theory, the chromatic polynomial is important for physics 
because of its connection with ground-state entropy. For 
a family G m , in the limit m —¥ oo (and hence n — > oo), 
we define 



W(Goo,<?) = lim P{G m ,q) 



l/n 



(1.4) 



2 



For a set of special q values, {q s }, the limits n — > oo and 
q — > q s do not, in general, commute for P(G, q) 1 ^, so one 
must specify the order of limits in defining W(G oa ,q) 
[H [Of. The set {q s } includes q = 0, 1, 2. For our 
purposes here, we define the order as q — > q s first and 
then n ~ > oo. For a large class of G m , if q is suffi- 
ciently large, then the number of proper g-colorings of 
G m grows exponentially with m and n, so that the Potts 
AFM has nonzero ground-state entropy per vertex on 
G oo ,S(G oo ,q) = k B ln[W(G oo ,q)]^0. 

The Sicrpinski graph S m is a planar graph constructed 
in the following iterative (i.e., hierachical) manner [![. 
One starts at stage m = with a triangle, then inserts 
vertices on the three edges and joins these with new edges 
to get Si. The central triangle is left unchanged, and 
at the next stage, m = 2, one repeats this procedure 
with the three other interior triangles, and so forth for 
higher to. The resultant graph at stage to is denoted S m . 
We note some basic properties of S m . The numbers of 
vertices and edges on S m are, respectively, 



n(S m ) 



3(3 r 
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(1.5) 



The degree of a vertex is defined as the number of edges 
that connect to it. The graph S m has three vertices of 
degree two, namely the ones at the three corners of the 
outermost triangle, while the other (3/2) (3 m — 1) vertices 
have degree 4. For a graph G with vertices of different 
degrees, it is convenient to define an effective or average 
vertex degree. This is 



K eff{G) 



2e{G) 



(1.6) 



Defining n e ff(G 00 ) = limn^oo 2e(G)/n(G), we have 

/W/OSoo) = 4 ■ (1-7) 

The chromatic number is easily shown to be x(Sm) = 3, 
and S rn is tripartite, i.e. 



P(S m ,3) = 3' 



(1.8) 



The girth of a graph G, g(G), is defined as the number 
of edges in a minimal-length closed (non-backtracking) 
path on G. This is g(S m ) = 3 for all m. The fractal 
dimension df of 6*00 is [l[ df(S OCl ) = In 3/ In 2 ~ 1.585. 
Some graphical properties of S m , including the numbers 
of spanning trees, spanning forests, dimer coverings, and 
Hamiltonian walks, have been studied recently in |14j |. 



II. CHROMATIC POLYNOMIAL AND ZEROS 

As with other hierarchical families of graphs, one can 
calculate Z(S m +i,q, v) from Z(S m , q, v) [9j, [l0|. We have 
done this up to m = 6. We first discuss the special case 
of P(S m ,q) = Z(S m ,q,—l). I n g eneral, if a graph G 
contains a complete graph K r [19[ as a subgraph, then 



lm(q) 
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FIG. 1: Zeros of the chromatic polynomial P(Si, q) in the q plane 
(123 zeros). 



P(G,q) contains the factor P(K r ,q) = IJ^o^ ~ j) 
Since S m contains one or more triangles K3 as sub 
graph(s), P(S m ,q) contains the factor q(q — l)(q — 2) 
From Eq. (fl~g]) . it follows that 



W(S 00 ,3) = 1 



(2.1) 



For q — 4, from our calculations up to m — 5, we infer 
that 



P(^ m ,4) = 2 3 ( 3m " +1 )-3 3m " 



(2.2) 



In the limit m -> 00, this gives W(Soo,^) = 2 2 / 3 • 3 2/9 = 
2.026346.., and hence the q — 4 Potts antiferromag- 
net has a nonzero ground state entropy [l(| S'(5 , 00 ,4) = 
(2/3)fc s [ln2 + (l/3)ln3] = (0. 706234.. )k B - 

The calculations of P(S m ,q) for m = and m = 1 
are elementary and yield P(So,q) = q(q — l)(q — 2) and 
P(Si,q) = q(q — l)(q — 2) 4 . We have computed the zeros 
of P(S m , q) in the q plane for m up to 5. We show plots 
of these chromatic zeros for m = 4 and to = 5 in Figs. 
[T] and [21 For to = and m = 1, the only chromatic 
zeros are real and occur at q — 0, 1, 2. For m = 1, the 
chromatic zero at q = 2 has multiplicity 4 and for the 
higher m values, we find that this zero has multiplicity 
3to. For all to values that we have studied, we find that 
the only real chromatic zeros of S m are at q — and 5=1 
(both with multiplicity 1), q = 2 with the multiplicities 
given above, and, for to > 2 at a value slightly less than 
q = 3, which monotonically increases toward q = 3 as m 
increases. Specifically, for to = 2,3,4,5, this largest real 
zero of P{S m ,q) occurs at 2.865331, 2.956063, 2.990362, 
2.9977175, respectively 

As in earlier work [17| , we denote the continuous accu- 
mulation set of chromatic zeros as to — > 00 as B q and the 
maximal point at which B q intersects the real axis as q c . 
As is evident in Figs. Q] and [51 the outermost zeros are 
consistent with the inference that in the to — > 00 limit a 
subset of zeros forms a closed loop on B q that crosses the 
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FIG. 2: Zeros of the chromatic polynomial P(Ss, q) in the q plane 
(366 zeros). 

real axis at q = and at a maximal point 

q c (Soo) = 3 . (2.3) 

For a one-parameter family of graphs G m , the locus B q is 
a slice, for constant v = — 1, of the locus in the C 2 space 
formed by the merging of zeros of Z(G m ,q,v) [l8j|. In 
general, if q c {G oo) G Z + , this corresponds to the property 
that for this value of q, B v passes through v = — 1, i.e., 
By passes through y = 0, signifying that the Potts AFM 
with q = q c {G oo) has a T = critical point on Goo. Eq. 
(|2.3[) is equivalent to the property that the q = 3 Potts 
AFM is critical on at T = 0. 

It is intriguing that, despite the fact that Soo shares 
three properties with the triangular lattice po| . namely 
(i) the girth equal to 3, (ii) the tripartite property, and 
(iii) the property of containing infinitely many triangle 
subgraphs, the value q c (Soo) = 3 is different from the 
value q c (tri) = 4 for the triangular lattice. Indeed, the 
value gc(S'oo) = 3 is the same as the value q c (sq) = 3 for 
the square lattice [HI . We showed via exact calculations 
that q c = 3 for infinite-length, finite- width self-dual strips 
of the square lattice [21] . The result qc(Soc) = 3 = q c (sq) 
is also interesting in view of the differences between Soo 
and the square lattice. Even before taking the m — > oo 
limit, these differences are evident. Thus, as noted above, 
S m has x{S m ) = 3 so that P(S m ,2) — and is tripar- 
tite, whereas a finite L x x L y section of the square lattice 
with free boundary conditions or toroidal boundary con- 
ditions with even L x and L y is bipartite, so for both of 
these families of graphs, denoted sq[L x x L y , BC], \ = 2 
and P{sq[L x x L y ,BC], 2) = 2. Furthermore, P(S m ,3) 
is a constant, namely 3!, while P(G[sq, L x x L y , BC},3) 
grows exponentially rapidly with n — L x L y . One prop- 
erty that Soo and the square lattice have in common is 
that K e ff(Soo) = K(sq) = 4. 

We remark further on the pattern of chromatic zeros 
in Figs. [T] and [21 The zeros on the upper and lower part 
of the outer loop exhibit a rather smooth appearance, 



although undulations on a sufficiently small scale would 
not be visible with our numbers of zeros. On the left side 
there is a decreased density of zeros, and on the right 
side the pattern of zeros has a more wavy and irregular 
structure. It is possible that some of the irregularly dis- 
tributed zeros on the right might form further continuous 
substructures on B q in the m — > oo limit. One also sees 
zeros in the interior that are consistent with the infer- 
ence that in the m — > oo limit they merge to form an 
inner closed loop on B q that crosses the real axis on the 
left at q ~ 2.62 and on the right at q ~ 2.74. We note 
that the left crossing point is equal, to within the uncer- 
tainty, with (3 + a/5)/2 = 4cos 2 (tt/5) = 2.6180.... There 
is also a substantial set of complex zeros scattered in an 
area whose center, on the real axis, occurs at q ~ 2.5; 
again, these could form further continuous substructures 
as m — > oo. 



III. GENERAL ZEROS OF Z(S m ,q,v) IN THE q 
PLANE 

As noted above, the chromatic polynomial is equal to 
the partition function of the Potts antiferromagnet at 
T = 0, i.e., v = — 1. It is of interest to study the zeros 
of the Potts partition function Z(S m , q, v) in the q plane 
for nonzero temperature in the antiferromagnetic interval 
v G (—1,0] and also for the ferromagnetic interval v > 0. 
We report these results here. For our purposes, it will 
suffice to use our calculation of Z(S^ q, v). 

In Fig. [3] we show zeros of Z(Si,q,v) in the q plane 
for the finite-temperature AFM value v = —0.5. These 
zeros are consistent with the inference that in the m — > oo 
limit, B q forms a closed locus that separates the q plane 
into two regions, crossing the real axis at q = and a q c 
value that has decreased from its value of 3 at v = — 1 
to ~ 1.5. (Strictly speaking, one cannot rule out the 
possibility of additional very small sliver regions, as were 
found in (22j, but we see no evidence for them here.) 
As was the case for regular infinite-length, finite-width 
lattice strips (e.g., [H, H1])j we nn d that q c decreases 
monotonically from its value of 3 at v = —1 to as v 
increases from —1 to 0. This is understandable, since for 
the Potts model on an arbitrary graph G, Z(G, q, 0) = q n , 
so that in the limit v — > 0, all of the zeros in q occur at the 
single value q = and the locus B q degenerates to this 
point. A notable structural change in these q-plane zeros 
that occurs as one increases v from —1 to —0.5 is that 
the zeros on the right side exhibit a smoother appearance, 
without visible outer or inner protruberances. 

We next consider the FM interval and show zeros of 
Z{Si, q, v) in q for the illustrative value v = 0.5 in Fig. S) 
These are consistent with the inference that for m — > oo, 
B q crosses the real axis at a value q ~ — 1.9 and does not 
cross this axis at any positive value. The latter result is 
in accord with the property that the q-state Potts FM on 
Sao has no finite-temperature critical point 0, 3 ■ Just 
as the density of zeros on the left side decreases in the 
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FIG. 3: Zeros of Z(Si,q,v) in the q plane for v = —0.5. 



FIG. 5: Zeros of Z(Sa, q, v) in the y plane for q = 3. 
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FIG. 4: Zeros of Z(S4,q,v) in the q plane for v = 0.5. 



vicinity of q — in the plots for v = — 1 and v = —0.5, 
so also for this FM value, v — 0.5, the density of zeros 
decreases in the vicinity of q = Q. 




FIG. 6: Zeros of Z(Sa, q, v) in the y plane for q = 5. 



IV. ZEROS OF Z(S m ,q,v) IN THE v PLANE 

We have also calculated zeros of Z(S m ,q, v) in the com- 
plex plane of a temperature-like Boltzmann variable v, 
or equivalently, y = v + 1 as a function of q. Plots of 
these zeros for q = 2, 3, 4, were presented in [l2j], and 
we confirm these results. Plots for larger values of q are 
also of interest, and we present these here. We denote 
the continuous locus formed by the accumulation set of 
these zeros in the limit m — > oo as B y . We show zeros 
of Z(S4, q, v) in the y plane for q = 3, 5, 10 2 , and 10 3 in 
Figs. 151181 

From our calculations of zeros of Z(S mi q,v) in v for 
various q and m, we can draw a number of inferences 
concerning B y for S^. First, for the Ising value q = 




FIG. 7: Zeros of Z(Si,q,v) in the y plane for q = 10 2 
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FIG. 8: Zeros of Z{Si, q, v) in the y plane for q = 10 3 . 

2, By passes through the point l/y = 0, in agreement 
with the fact that the Ising FM has a T — critical 
point on the Sierpinski gasket Soo , as evidenced from the 
divergence in the correlation length as T — > • In the 
l/y plane, Bi/ y forms roughly a figure-8 shape, crossing 
the imaginary l/y axis at ±0.7i. Second, equivalent to 
our conclusion from the chromatic zeros that q c (S oo) = 3, 
we infer that for the q — 3 Potts AFM on , B y passes 
through y = 0. Our results for q c (Soa) also imply that 
By does not pass through y = for the g-state Potts 
antiferromagnet with q > 3. Our calculations of zeros 
in y for q > 3 agree with this. As q increases beyond 
q = 3, the inner loop of zeros moves to the left. As 
is evident in Fig. [8j for sufficiently large q, there is no 



longer any visible inner loop. As q gets large, the outer 
loop of zeros becomes more circular, and we infer that 
in the limit m —> oo this is also true for B y . Since the 
o-state Potts ferromagnet is critical at T = on 
Q, it follows that B y always passes through the point 
l/y = 0. As is evident from Fig. this means that as one 
approaches the real axis on B y , the locus deviates from 
the roughly circular pattern to pass through l/y = 0. 
These inferences are in accord with our earlier work in 
[23 | (see also [2^], where it was shown that the radius of 
this approximately circular pattern of zeros is 

\y\ = q 2/Ke " as q -> 00 . (4.1) 

Since K e // approaches 4 on S m with increasing m, this 
radius approaches \y\ = Jq- Hence, this radius is approx- 
imately 10 for the case q = 10 2 in Fig. [7j and ~ 32 for 
q = 10 3 in Fig. [5J in agreement with our results shown 
in these figures. 

We have also carried out similar calculations for gen- 
eralizations of the Sierpinski gasket, including the gen- 
eralization to three dimensions. With R. Roeder, we 
have performed an analogous study for another family 
of graphs, the m'th iterates of the diamond hierarchical 
lattice, D m , which lead, in the m — > 00 limit, to the frac- 
tal Dqo. We find that for D m , the chromatic zeros are 
considerably more distributed throughout regions of the 
q- plane than is the case with S m . Our results will be 
reported elsewhere. 
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